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Abstract: We show how to construct consistent braneworld models which exhibit 
late time acceleration. Unlike self- acceleration, which has a de Sitter vacuum state, 
our models have the standard Minkowski vacuum and accelerate only in the presence 
of matter, which we dub "stea/t/i-acceleration". We use an effective action for the 
brane which includes an induced gravity term, and allow for an asymmetric set-up. 
We study the linear stability of fiat brane vacua and find the regions of parameter 
space where the set-up is stable. The 4-dimensional graviton is only quasi-localised in 
this set-up and as a result gravity is modified at late times. One of the two regions is 
strongly coupled and the scalar mode is eaten up by an extra symmetry that arises in 
this limit. Having filtered the well-defined theories we then focus on their cosmology. 
When the graviton is quasi-localised we find two main examples of acceleration. In 
each case, we provide an illustrative model and compare it to ACDM. 
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1. Introduction 

Recent observations of high red-shift supernovae suggest that our universe has en- 
tered a period of acceleration [T[|. Combining these observations with the microwave 
background data, 0, and the results of large scale structure surveys P], we deduce 
that the universe is approximately spatially fiat, and that the observed acceleration is 
consistent with standard AD cosmology if we assume that roughly 70% of the actual 
energy content of the universe comes from a small positive cosmological constant, 
A/SttG ~ 10~^^ (eV)'^ 0]. Of the residual matter content, only about 4% corre- 
sponds to baryonic matter, the remaining component corresponding to an also as yet 
unobserved "dark matter sector" . Although this simple energy pie chart fits the data 
extremely well, effective field theory methods have so far failed to give a satisfactory 
explanation for such a small value of the cosmological constant. Naive arguments 
predict a cosmological constant of the order of mpj, i.e. 10^^° times larger than the 
"observed" value, which fine tuning problem has inspired a search for alternative 
explanations for the observed acceleration. 

An interesting alternative is that supernovae data are actually indicating the 
presence of a new physical scale where novel gravitational physics kicks in. The scale 
in question is a classical energy scale corresponding to the current Hubble curvature 
radius, Hq ~ 10^'^^ (eV). The braneworld scenario ^ has been a natural breeding 
ground for these ideas (see, for example, [^, ^, [13 )• Typically, gravity becomes 
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weaker in the far infra-red due to gravitational leakage into the extra dimensions. In 
some cases, this modification of gravity at large distances leads directly to exponential 
acceleration at late times, even when there is no effective cosmological constant P, 



T0| , |TT[. These types of solutions are often referred to as self- accelerating, since 
their vacuum state corresponds to a de Sitter brane. However, there is by now 
good evidence that self-accelerating solutions lead to problems with perturbative 
ghosts [0. 

Modifying Einstein's General Relativity in a consistent manner in the infra- 
red has proven to be a very difficult task. There are three main reasons for this: 
theoretical consistency such as stability or strong coupling; experimental bounds 
mostly coming from the solar system; and naturalness, i.e. a restricted number of 
additional fields and parameters for the modified theory. The first problem, for the 
time-being at least, has been the most stringent one. Indeed, it is common for the 
vacua of such modified theories to suffer from classical and/or quantum instabilities, 
at least at the level of linear perturbations. Most candidate braneworld models, 
[0, ^, ^ have been shown to suffer from such instabilities or strong coupling or both, 
0, |l3l. Generically, a ghost mode appears in the perturbative spectrum of the 
theory at the scale where gravity is modified, effectively driving the acceleration, 
much like a cosmological constant put in by hand. It has been argued that strong 
coupling, i.e. the breakdown of linear perturbation theory |15[, may invalidate 
these claims, bypassing also rather stringent solar system constraints. As a result 
however, in the absence of exact solutions the gravity theory becomes unpredictable 
at the strong coupling scale and the loss of linear perturbation theory is perhaps a 
bigger problem to overcome. 

The situation is not better settled for four dimensional theories that modify the 
Einstein-Hilbert action directly. Here, early attempts go back to Brans and Dicke 
who considered a massless scalar-tensor theory with only one additional con- 
stant degree of freedom, the kinetic coupling ubd- Such a theory (which for example 
goes through the naturality condition) has been ruled out by solar system time delay 



experiments pushing ubd ^ 40000 coming from the Cassini satellite |T^. Non- 



minimal scalar-tensor theories with some potential and hence varying u can still be 
consistent with PPN experiments if we suppose that the scalar field varies between 
solar system and cosmological scales. Indeed, we remind the reader that the rele- 
vant cosmological scale Hq^ and the size of the solar system give a dimensionless 
quotient of {IAUHq) ~ 10~^^ whereas the typical PPN parameter 7 is 1 with error 
10~^. Vector-tensor theories have also been shown to be plagued by instabilities (the 
exception being the case studied by Jacobson |18[ where the vector has a fixed norm, 
and Lorentz-invariance is broken at the level of the action^). Adding higher order 
combinations of the Riemann curvature tensor also generically leads to the appear- 

-•^At the level of cosmological solutions this is of course true since there is no timelike Killing 
vector; the difference there being that the symmetry is not broken at the level of the action. 
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ance of ghosts except if one considers functionals of the Ricci scalar curvature F{R) 



W^ . Then the theory in question can be translated via a conformal transformation 



to a non-minimal scalar-tensor theory. Solar system constraints are again problem- 
atic (though see ||2^ for recent developments) and often vacuum spacetimes such as 
Minkowski are not even solutions to the field equations (although perhaps this is not 
such a big problem in the cosmological setting). 

In this paper we present an alternative general class of ghost-free braneworld 
models which possess accelerating solutions. However, unlike the usual self-accelerating 
braneworld solutions, and indeed unlike ACDM, our models do not necessarily lead 
to de Sitter cosmologies at late times. Indeed, the final fate of the universe may 
not even be an accelerating cosmology! Also, unlike the usual models in which even 
the vacuum brane accelerates, in our models the vacuum state is a stable Minkowski 
brane, and the Friedmann equation actually implies that while the Minkowksi brane 
is the vacuum state for an empty brane, once one has an "ordinary" cosmological 
fluid, i.e. p = wp with w > 0, then the brane can enter an era of acceleration, even 
though there is no cosmological constant or dark energy present. Depending on the 
specific parameters in the model, this acceleration can either persist, or the universe 
at some stage exits from accelerated expansion into a stiff-matter cosmology. In 
either case, the expansion is power law, rather than exponential. 

Clearly an accelerating cosmology is of no interest if it is hampered by the usual 
consistency problems, therefore we perform a general scan of possible braneworld 
solutions, filtering out those solutions that we suspect will contain problems with 
ghosts (such as self-accelerating vacua or regions where we prove the appearence of 
ghosts). We also further restrict our parameter space by requiring that the bout 
of acceleration is preceded by a sufficient period of standard AD cosmology, in or- 
der to reproduce standard early universe cosmology, such as nucleosynthesis, the 
development of structure etc., so that the standard cosmological picture is retained. 
Note therefore, that these braneworld models are not suited for explaining primor- 
dial inflation. In fact, the basic feature of these models is that it is only at late 
times in a matter dominated universe that this acceleration can kick in. Indeed, if 
the parameters of the model were tuned so that the modified gravity scale were to 
become important during the radiation era, then the universe would not accelerate, 
but merely 'coast' [^, before returning to an expanding, decelerating, cosmology. 
Therefore, our models explain quite naturally why it is only at scales comparable to 
the current Hubble time (or indeed longer) that acceleration can occur. 

The rest of this paper is organised as follows. In the next section we describe our 
set-up involving a single brane embedded in some five- dimensional bulk spacetime. 
This will be a combination of the DGP 0], and asymmetric brane [10| models. In 
other words, we consider an asjTumetric set-up, with an induced curvature term on 
the brane. We will present the equations of motion, and derive the vacuum states 
corresponding to Minkowski branes. In section ^ we consider vacuum perturbations. 
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focusing on the scalar radion mode. We will derive the effective action for this mode, 
in order to filter out those solutions that contain a perturbative ghost. In section ^ 
we introduce a cosmological fluid and derive the cosmological solutions, establishing 
the conditions required for a consistent cosmology evolving towards a Minkowski 
vacuum as the energy density decreases. The consistent solutions are then analysed 
in greater detail in section ^, with special emphasis on the conditions for cosmic 
acceleration. We analyse two main models in detail, showing how they give rise to 
cosmic acceleration, and comparing these models with ACDM. Finally, we conclude 
in section ||. 



2. Asymmetric braneworlds: formalism and set-up 

Consider a single 3-brane, S, embedded in between two bulk five-dimensional space- 
times Aii, where i = 1,2. The brane can be thought of as the common boundary, 
S = dAii = dM.2 of these manifolds. Each spacetime M.i generically has a five 
dimensional Planck scale given by Mj, and a negative (or zero) cosmological con- 
stant given by Aj = —Qkf. In general we will not consider Z2 symmetry and we will 
allow for the cosmological constants, and even the fundamental mass scales to differ 



on either side of the brane JTOl- Allowing for Ai 7^ A2 is familar enough in domain 
wall scenarios Here we are also allowing for Mi 7^ M2. This is not so familiar, 
but could arise in a number of ways. Suppose, for example, that this scenario is 
derived from a fundamental higher dimensional theory. This theory could contain 
a dilaton field that is stabilised in different fundamental vacua on either side of S. 
From the point of view of a bD effective description, the 5D Planck scales would then 
differ accordingly. Indeed naive expectations from string theory point towards this 
asymmetric scenario as opposed to a symmetric one. Different effective Planck scales 
can also appear on either side of a domain wall that is bound to a five-dimensional 
braneworld |2^. 

The brane itself has some vacuum energy, or tension. This will ultimately be 
fine-tuned against the bulk parameters, Aj and Mj, in order to admit a Minkowski 
vacuum solution. In other words, there will be no effective cosmological constant on 
the brane. As in the original DGP model we will also allow for some intrinsic 
curvature to be induced on the brane. Such terms are rather natural and can be 
induced by matter loop corrections p3, finite width effects or even classically 



from higher dimensional modifications of General Relativity 

Our set-up is therefore described by the general 5 dimensional action, 

S = •S'bulk + 'S'brane (2-1) 

The bulk contribution to the action is described by 

5buik =y^Mf f v^(i? - 2A,) + 2M^ [ y/^K^'^ (2.2) 
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where gab is the bulk metric with corresponding Ricci tensor, R. The metric induced 
on the brane (dAii) is given by 



Jab = 9ab - riaUb (2.3) 

where is the unit normal to d}Ai in M.i pointing out of AAi. Of course, continuity 
of the metric at the brane requires that 7ab is the same, whether it is calculated from 
the left, or from the right of the brane. In contrast, the extrinsic curvature of the 
brane can jump from right to left. In dM.i, it is defined as 

K^i = llli^icu,) (2.4) 



Its trace appears in the Gibbons- Hawking boundary term in 

The brane contribution to the action, meanwhile, is described by 

>S'brane = / V-^i^llT^ - CT + ^matter) (2.5) 
J brane 

where a is the brane tension, and /^matter includes any matter excitations. We have 
also included the induced intrinsic curvature term, 71, weighted by a AD mass scale, 
m-pi. Note that we have taken > 0, as in the original DGP model. There are two 
reasons for this. Firstly, when it comes to studying the cosmological solutions, we 
would like this term to dominate the cosmology at early times, in order to reproduce 
the standard 4Z) cosmology, as discussed in the introduction. Secondly, allowing 



m 



pi 



< could result in vacuum perturbations containing a spin-2 ghost [|2^ . 



The equations of motion in the bulk region, Aii, are just the Einstein equations, 
with the appropriate cosmological constant, Aj. 

Eab = Rab - ^Rgab + ^i9ab = (2.6) 

The equations of motion on the brane are described by the Israel junction conditions, 
and can be obtained by varying the action (|2.1|) , with respect to the brane metric, 
-fab- This gives 

eab = 2 {M^iKab - K-fab)) + mj, (^ab - ^7^7a6^ + '^lab = ^Tab (2.7) 

where Tab = — ^^^- ^^t"""" • The angled brackets denote an averaged quantity at 
the brane. More precisely, for some quantity Qi defined on the brane in dAii, we 
define the average 

(Q) = . (2.8) 

Later on we will also make use of the difference, AQ = Qi — Q2- Note that the 
Israel equations here do not use the familiar "difference", because we have defined 
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the unit normal as pointing out of A4i on each side, i.e. the approach is that of the 
brane as a boundary. Israel's equations on the other hand were derived for thin shells 
in GR, i.e. where the brane is a physical, very thin, object, and the normal is thus 
continuous, pointing 'out' on one side of the wall, and 'in' on the other. 

We will now derive the vacuum solutions to the equations of motion (|2.6|) and 
( p.7|) . This corresponds to the case where there are no matter excitations, and so. 
Tab = 0. In each region of the bulk, we introduce coordinates x"" = {x^,y), with the 
brane located at ?/ = 0. We will not be interested in de Sitter solutions, since these 
will only arise through an excess in vacuum energy, or through "self-acceleration". 
The former offers no alternative to ACDM, whereas the latter is expected to suffer 
from a generic ghost-like instability [|1^]. Therefore, we seek solutions of the form 



ds"^ = gabdx°'dx^ = a^{y)rjabdx"-dx'' (2.9) 
Inserting this into the bulk equations of motion ( p. 61) gives 



-V = fcV, - = 2A;V (2.10) 
a J a 

where "prime" denotes differentiation with respect to y. Note that we have dropped 
the asymmetry index i for brevity. Equations ( |2.10 ) have solution 



where 6* = ±1. Note that each region of the bulk corresponds to < ?/ < ^max, where 

.Ilk for^ = l , ^ 

ym.. = { ' (2.12) 
oo for 6* = —1 



For k 0, this means that when ^ = 1 we are keeping the adS boundary (growing 
warp factor) whereas when 6' = — 1 we are keeping the adS horizon (decaying warp 
factor). For A; = 0, we simply have a Minkowski bulk, in the usual coordinates, and 
the sign of 6 is irrelevant. 

The boundary conditions at the brane (|2.7] ) lead to a finely tuned brane tension 



a = -12{M^ek) (2.13) 

This fine tuning guarantees that there is no cosmological constant on the brane, and 
is equivalent to the (asjTiimetric) Randall- Sundrum fine-tuning for 6iki < and 

e2k2 < §. 
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3. Linearised vacuum perturbations and asymptotic stability 



We shall now consider linearised perturbations, hab about our background solutions 
Q and ( PH ), so that 

ds^ = a\y) [r]ab + hab{x,y)] dx^dx^ (3.1) 

In the unperturbed spacetime, the gauge was fixed in both Jvli and M.2 so that the 
brane was at y = 0. However, a general perturbation of the system must also allow 
the brane position to fiutter. In A^j, the brane will be located at 

y = Mxn (3.2) 

Of course, these expressions contain some gauge-dependence due to invariance under 
the following diffeomorphism transformations 

y^y + r]{x,y), x^' ^ x>' + C^{x,y) (3.3) 

Now, it is convenient and physically relevant to decompose these transformations in 
terms of the AD diffeomorphism group. This gives = ^'^ + d'^C,, where is a 
Lorentz-gauge vector satisfying c?^^'^ = 0. We do likewise for the perturbation: 

V = + 29(^F,) + 2d^d,E + 2Ar]f,, (3.4) 
hf,y = B^ + d^B (3.5) 
hyy = 2(f) (3.6) 

Again, and are Lorentz-gauge vectors, whereas /i^J is a transverse-tracefree 
tensor, d^h^^^u = h^'^^^^ = 0. Note that greek indices are raised and lowered 
using rjfj^^. Under the gauge transformations (|3.3|), the various components of the 
perturbation transform as follows 

hll - (3.7) 
B^^ B^- F^- (3.8) 

B^B-C-V, E-^E-^, A-^A--r] (3.9) 
a a 

We immediately see that the tensor component, /i^J is gauge invariant. We can also 

construct the following vector and scalar gauge-invariants in the bulk 

X, = B,-F'^, X = A-^{B-E% r = 0-M^^i:^ (3.10) 

We will now consider vacuum fluctuations, so that the bulk equations of motion and 
the Israel junction condition are given by 6Eab = SQab = 0. At this point we assume 
that the tensors, vectors and scalars do not mix with one another, so that their 
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equations of motion can be taken independently. We will come back and check the 
validity of this hypothesis against our results later on. Focusing on the scalars, we 
find that the bulk equations of motion give SEH^^^"^ = 0, where 



^^scalar 



dy + 3- 



X' 



-Y 



6E: 



iscalar 
MS/ 



X' Y 

a 



6E, 



scalar 

yy 



Sd^X + 12- X' 



-Y 



(3.11) 



(3.12) 



(3.13) 



These equations are easily solved on each side of the brane (we remind the reader 
that we have dropped the asymmetry index i) to give 



X 



U{x) 



Y 



2U{x) 



(3.14) 



where d'^U = 0. We can identify U{x) as the bulk radion mode. It represents two 
degrees of freedom, Ui{x), i = 1,2, one on each side of the brane. In addition to this, 
we have the brane bending degrees of freedom, given by /i and /2. Given that we 
have two boundary conditions at the brane, namely continuity of the metric and the 
Israel equations, we expect the number of physical scalar degrees of freedom to be 
at most two, as we shall now demonstrate. 

To impose the boundary conditions at the brane it is convenient to work in 
brane-GN gauge. In this gauge, we have B = (p = and the brane is fixed at y = 0. 
It follows that 



E = W{x) + V{x)Q{y)-U{x)Q{yf 



A 



U{x) 



[Vix)-2U{x)Qiy)] (3.15) 



where 



Q{y) 



dyi 



[1 



-2 



(3.16) 



The brane bending degrees of freedom are now encoded in the fields, Vi and V2. In 
contrast, Wi and W2 merely reflect the freedom to choose the gauge along the brane, 
and as such, can be consistently set to zero. To see this, we evaluate our solution at 
y = to derive the scalar part of the brane metric ( |3.4| ) 



scalar , 



X, 0) = 2d^d,W + 2[U - ekV]ri^, 



(3.17) 



The pure gauge part, 2d^dyW , and the remainder 2\U — 6kV]rj^i, need to be well 
defined independently of one another. This means that Aiy = 0, and it is easy to 
to see that we can continously gauge away W on both sides of the brane, so that 
W = Q. 
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The continuity condition at the brane and the Israel equations now require that 



A[U - 9kV] = {M^V + mli[U - OW]) = (3.18) 

Therefore, as predicted, we are left with at most two scalar degrees of freedom as 
perceived by a 4 dimensional observer. Actually, it is easy to see why this is so if 
one introduces regulator branes at (say) y = y^, < 2/max- The remaining degrees of 
freedom essentially correspond to the fluctuations in the proper distance in between 
branes. Actually, the introduction of regulator branes is particularly useful when 
trying to impose boundary conditions in the bulk, and in terms of calculating a AD 
effective action. We shall now develop this in more detail. 

In the background, the extrinsic curvature, K*^^, of the regulator brane satisfies 
the following 

a?{y^ 

where 7*^ is the induced metric on the regulator brane. We require that this equation 
also holds in the perturbed scenario. This time, it is easiest to work in GN coordinates 
relative to the regulator brane (we shall call these regulator-GN coordinates). The 
regulator brane is now fixed at y = y^,, and B = (p = 0. Again, it follows that 

E = W,{x) + V,{x)Q,{y) - U{x)QM\ ^ = - 4 [^*(^) - '^U{x)QM] 

(3.20) 

where 



As before we can set = 0. Applying the boundary condition (|3.19|) we deduce 
that K = 0. 

It is very important to realise that the brane-GN coordinates and the regulator- 
GN coordinates are not necessarily the same. In particular, in brane-GN coordinates 
the true brane is at ?/ = 0, but the regulator brane might not he at y = y^. Of course 
the two sets of coordinates must be related by a gauge transformation. Starting in 
regulator-GN gauge, we transform to brane-GN gauge with the following coordinate 
change 

x^ -^x'' + D^C, y ^y + v (3.22) 

where 

e = U [Q{yy - QM'] - VQ{y), ^ = 7 - ^ ^^^^^ ~ ^^^^^^ ^^'^^^ 

From this we can deduce that in brane-GN coordinates the regulator brane is at 
y = y* + vi^-, y*) where 

(3.24) 

ayy*) 
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We would now like to calculate the proper distance between the branes. To do this 
carefully, it is best to work in a gauge in which both the true brane and the regulator 
brane are fixed. Clearly neither brane-GN nor regulator-GN gauge will do. Instead, 
we start off in brane-GN coordinates, and transform to fixed brane coordinates by 
letting 

x^^x'', y^y + e{x,y), (3.25) 



so that 



(ea)' „ „ „ . , a' 



B ^ -e, E-^E, A^A e (3.26) 

a a 

If we set e(x, 0) = and e(x, y*) = —rj{x, y*), both branes will be fixed, at ?/ = and 

y = y* respectively. The proper distance between them is given by 

ry* 

z* = dy^/g^i^ = + 6z^ (3.27) 
Jo 

where ^ 

= -—\n{l - Oky,) (3.28) 

is the proper distance between branes in the background, and 

6z, = V-2Q{y,)U (3.29) 

is the fiuctuation. Fixed wall gauge is also useful for calcinating the effective action. If 
we assume that the boundary conditions hold at the true brane and at both regulator 
branes, the effective action (to quadratic order) is simply given by the following bulk 
integral 

Ses = I /b,ik V^M^g-'SEa, (3.30) 

Choosing fixed brane gauge in both the left hand bulk and the right hand bulk, we 
find 

S^s= j d^x(M'^5z,d'^U{x)^ (3.31) 

where U can be related to the fields 5z^ using the boundary conditions at the brane 
( [3.181 ). Now consider what happens, as we remove the regulators by taking the limit 



y 1/max, or equivalent ly — > cx). For U ^ 0, 

{—A/k whenever 9k > , , 

' (3.32) 
oo whenever 6k < 

Therefore, if 6'A; < 0, the fiuctuation in the proper distance between branes diverges. 
This refiects the fact that the gauge invariant mode, U, is non-normalisable. Since we 
must have normalisable boundary conditions in order to obtain a local 4D effective 
theory on the brane, we require that U = whenever 6k < 0. The result is that 
the scalar degree of freedom coming from a given bulk region only survives the 



- 10 - 



single brane limit when the corresponding regulator is taken to the AdS boundary as 
opposed to the AdS horizon. This makes sense for the following reason. If an observer 
on the brane shines a light ray into a bulk region towards the AdS boundary, he/she 
sees it reflected back from the boundary in a finite proper time. However, if the ray 
is shone into the bulk towards the AdS horizon, it will never be reflected back. One 
is therefore able to detect fluctuations in the proper distance between the brane and 
the AdS boundary, but not the brane and the AdS horizon. 

After removing the regulators, and imposing the boundary conditions (|3.18| ), we 
find that there are no normalisable scalar degrees of freedom left if 6iki,62k2 < 0. 
This does not come as a surprise since in this case the bulk includes the AdS horizon 
on both sides of the brane, and so one would not expect there to be any normalisable 
radion degrees of freedom for the reasons outlined in the previous paragraph. In all 
other cases, the effective action ( p. 311 ) is given by 



Sea = 6 



d^x{dAf (3.33) 



where the field A = A{x, 0) = U — dkV measures conformal rescalings of the metric 
on the brane. It represents the only scalar degree of freedom that survives, since it is 
the only scalar in the induced metric that cannot be locally gauged away. Physically, 
however, we might have expected there to be two radion degrees of freedom since both 
sides of the bulk included the AdS boundary. Careful inspection however, shows that 
this is not the case. Indeed reintroducing the regulator branes, we can split the two 
degrees of freedom into a centre of mass motion of the brane and a relative motion 



2S]. The latter is shown to drop out as we push the regulators out to the boundary. 



Furthermore, note that even the centre of mass motion drops out if we have a zero 
cosmological constant on either side of the brane. This is because in flat spacetime 
we have an extra translation Killing vector in our coordinate chart (with respect to 
adS) that permits us to consistently gauge away the remaining scalar mode. 

Before we start to derive conclusions about when there is a scalar ghost in the 
spectrum of perturbations, let us pause for a moment to discuss the reliability of 
these results. In deriving the scalar effective action we assumed that the scalars 
did not mix with the vectors or tensors. Now, we need not worry about vectors, 
because the vector contribution can always be locally gauged away on the brane. 
The same cannot be said for the tensors, which are gauge invariant, so there is a 
danger that they could mix with the scalars, as was the case for the self-accelerating 
DGP solution [|r^. To see whether this can happen here, we write the scalar gauge- 
invariant piece in Gaussian-Normal gauge as follows 

= -^A9^U (3.34) 

Because d'^U = 0, we can view massless transverse-tracefree tensor. There- 

fore, we can trust our analysis provided there are no massless modes in the tensor 



- 11 - 



sector. This is guaranteed if the background bulk has infinite volume. The volume of 
the background is finite if and only if 6iki and ^2^2 are both strictly negative. This 
corresponds to a generalised Randall- Sundrum scenario where the warp factor decays 
into the bulk. Then indeed there is a tensor zero mode, and we have mixing between 
tensors and scalars at zero mass. Although our analysis may therefore be unreliable 
in this instance, it is well known that in this case there is actually no radion degree 
of freedom. In fact, the mixing actually ensures that the graviton propagator has 
the correct 4D tensor structure |^ as in the original RS scenario. It also follows 
that there can be no infra-red modification of gravity since the graviton zero-mode 
always guarantees Einstein like behaviour at large distances and at late times. This 
will not be so interesting from an infra-red cosmological perspective, so from now on, 
we will drop the case 9iki, 6*2^2 < 0. 

In all other cases, the bulk volume is infinite and there is no normalisable zero- 
mode graviton in the spectrum and hence no mixing. However, we do want a quasi- 
localised graviton dominating up to some energy scale in the graviton spectrum as 
one gets for DGP |p or GRS 0. This is guaranteed by the induced curvature on 
the brane entering with the "correct" sign (mpj > 0), as we will demonstrate in 
the next section. Whenever the quasi-localised zero mode dominates, the higher 
mass modes and the radion mix giving some type of well-defined generalised scalar- 
tensor gravity. At large distances, however, the quasilocalised nature of the graviton 
disappears giving a continuum of massive modes with no mass gap. Then the radion 
mode ( |3.34| ) no longer mixes and is dangerous in the sense that it can be a ghost. 
Whenever kik2 = we have already discussed how the radion mode disappears. 
However, if kik2 7^ 0, the radion is present, and we see from ( p.33| ) that it is not a 
ghost, provided 



It is interesting to note the competition between the bulk term in x, and the induced 
gravity term: the "correct" sign for the DGP term (mpj > 0) always contributes to a 
ghost-like radion (as does a localised warp factor). This behaviour is caused by the so- 
called conformal ghost that appears in 4D Einstein gravity for perturbations about 
Minkowksi space (see for example |]30[). In that case the scalar mode is harmless 
since it provides the correct tensor structure in the propagator mixing with the zero- 
tensor mode. Strictly speaking each individual mode doesn't exist independently: it 
is only their linear combination which has physical meaning. However, in our case 
the absence of a graviton zero-mode ensures that there is no mixing between modes 
and so a radion ghost always signals a vacuum instability. 

Of particular interest are the limits x 0, and x ~^ 00. For small/large the 
radion couples very strongly /weakly to the trace of the energy- momentum tensor, 
since schematically we have. 




(3.35) 



Xd' 



(3.36) 
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In the X limit, the radion completely decouples. This corresponds to the case 

where kik2 = so at least one side of the bulk is Minkowski. The radion decouples 
in this case because it costs no effort to translate the brane toward the Minkowski 
side, and hence no brane bending can be detected. 

The X limit might also be refered to as the conformal limit, since the brane 
can only support conformal matter sources (T = 0) to linear order in perturbation 
theory. If we introduce some non-conformal matter, the linearised theory breaks 
down because of strong coupling, and the geometry responds non-linearly. This 
behaviour reflects the onset of a new symmetry. It is reminiscent of the partially 
massless limit (m^ = 2H^) of a massive graviton propagating in de Sitter space . 
In that theory an extra symmetry kicks in that eliminates the scalar degree of free- 
dom. In our case, the linearised field equations, SEab = and 6Qab = |^a6, become 

ab 't'ab "I" ft'ab ' ^ 



1 1\ 

invariant under the transformation hob — > hab + h„J where 



= (1 - a-')d,d^f - 2P/r],., h^f^ = hlfj = (3.37) 

This is pure gauge in the bulk, but not on the brane. The transformation therefore 
encodes an extra symmetry, beyond the usual diffeomorphisms, that eliminates the 
radion degree of freedom when x = 0- 

In summary, if the bulk has infinite volume, we have a theory without a nor- 
malisable tensor zero-mode, for which the radion will dominate at large distances. 
In certain limits, the radion either decouples (x —>■ oo), or is eliminated by a new 
symmetry (x — ^ 0). Otherwise, the radion is present, and will render the vacuum 
unstable on large scales unless x ^ 0. 



4. Cosmological solutions 

Let us now consider what happens when we introduce a cosmological fluid to the 
Minkowski braneworlds derived in section ^. In order to preserve homogeneity and 
isotropy on the brane, we must assume that the bulk metric is a warped product of 



the form 32 



ds^ = A\t, z){-de + dz^) + B\t, z)d^l (4.1) 

where for n = 1, 0, —1, rfx^ is the metric on the unit 3-sphere, plane, and hyperboloid 
respectively. It turns out that we have enough symmetry to render the bulk equations 



of motion (|2.6| ) integrable, and a generalised form of Birkhoff's theorem applies [p^j . 
When A = K = 0, the bulk solution is just a portion of Minkowski spacetime, whereas 
in all other cases the bulk is a portion of a black hole spacetime with cosmological 
constant, A = — 6fc^, and horizon geometry parametrised by k: 



ds^ = -V{r)dt^ + 7^ + r'd^l V{r) = kh-^ + k-^ (4.2) 
V[r) 
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In order to construct the brane, we glue a solution in Aii to a solution in A^2, 
with the brane forming the common boundary. Let us describe this in more detail. 
Assume for the moment that the bulk solution on both sides of the brane takes the 
form ( |4.2| ). In A4i, the boundary, dJ^i, is given by the section ()f:j(r), rj(r), x^) of 
the bulk metric. The parameter r is the proper time of an observer comoving with 
the boundary, so that 

where overdot corresponds to differentiation with respect to r. The outward pointing 
unit normal to dAii is now given by 

n, = e^iniT),-U{T),0) (4.4) 

where &{ = ±1. For di = 1, Aii corresponds to rj(r) < r < oo, whereas for 6i = — 1, 
Aii corresponds to ^ r < rj(r). The signs of 6 are consistent with the analysis in 
the previous two sections. 

The induced metric on dAii is that of a FRW universe, 

ds^ = -dr^ + Tiirfd^l, (4.5) 

Since the brane coincides with both boundaries, the metric on the brane is only well 
defined when ri(r) = r2(r) = R{t). The Hubble parameter on the brane is now 
defined hy H = R/R where R = R{t) is the brane trajectory in the bulk spacetime. 



The brane equations of motion are again given by the Israel equations (|2.7| ). In 
addition to the finely tuned tension, we introduce a homogeneous and isotropic fluid 
on the brane, with energy density, p, and pressure, p. We will assume that these 
satisfy the strong energy condition, so that p > and p + 3p > 0. This guarantees 
that any cosmic acceleration that may occur is entirely due to modified gravity. The 
Israel equations now give the following 



"^'-m%(H^ + ^)-2 ImHJh^ + ^ ) (4.6) 



where the tension, a, is given by equation ( p.l3| ). Note this formula is general and 



holds for all values of parameters, even those involving Minkowski backgrounds, for 
which we simply set V{R) = in ( [4.6| ). 

For simplicity, we will assume from now on that the bulk spacetime is maximally 
symmetric by setting the mass term, p = 0. One can confirm that this term, which 
has the behaviour of dark radiation, does not play an important role at late times, 
when it is subdominant to the matter content on the brane. In addition, since obser- 
vations of the first acoustic peak in the Cosmic Microwave Background demonstrate 
that the universe is very nearly flat 0, we set k = 0. The modified Friedmann 
equation then takes the compact form 

p = F {H') (4.7) 
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where 

F{H^) =QmliH^ -12 (^M^elyiPT^ -k^^ (4.8) 

Now, as if , it is easy to check that the energy density p —>■ 0. This serves as a 
consistency check that the finely tuned tension (p.l3|) guarantees that the Minkowski 
brane corresponds to a possible vacuum solution. There might be other vacuum 
solutions with Hq > 0, corresponding to self-accelerating vacua with F{Hq) = 0. 
These however are suspected to contain perturbative ghosts |jl2| (although we will 
not attempt to show this explicitly here). We will return briefiy to discuss how 
and when such vacua arise presently. Also note that switching off Mj we get the 
usual Friedmann equation, and switching off m^; gives the usual Randall- Sundrum 
modified Friedmann equation p3| . 

Since we wish to study the possibility of cosmic acceleration in these models, 
it is natural to examine the cosmic deceleration parameter, q = Assuming a 

constant equation of state, p = wp, it follows from (|4.7| ) and conservation of energy 

p = -3H{p + p) (4.9) 

that the deceleration parameter is given by 

q = -1 + ^(1 + w)C{H^) (4.10) 

Here the functional 

will vary during the course of the cosmological evolution whereas for GR we have 
simply C = 1. It follows that the deceleration parameter will also vary during the 
cosmological evolution, in contrast to what happens in four dimensional Einstein 
gravity where it is constant for constant w. For cosmic acceleration to occur, ob- 
viously the deceleration parameter must become negative. For ordinary forms of 
matter satisfying the strong energy condition 1 -|- 3w > 0, it is easy to check that 
acceleration can only be achieved if C{H'^) falls below one. 

It is useful to define an effective equation of state parameter for the universe, 
7eff, which is related to the deceleration parameter q = |(1 + 37efT), from which we 
deduce that 

l+7eff = C(l + w;) (4.12) 

Again, as C varies during the cosmological evolution, so must the effective equation of 
state. Note however, this should not be confused with a dark energy component with 



a varying equation of state (see e.g. |34]), in which there are two main components of 
the energy of the universe: matter and the time- varying dark energy. Here, we assume 
that there is only one dominant cosmological fluid, and 7eff provides a simple, effective 
way of tracking the gravitational effect of that fluid. It is clear that the cosmological 
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behaviour is completely determined by the functional, F, and its derivatives, with 
the combination C = F/ F'H^ proving particularly important when asking questions 
about cosmic acceleration. To emphasize the difference between taking dark energy, 
in which we modify the matter in the Friedmann equation to = J^{p), and light 
gravity in which we modify the geometry: p = F^H"^), consider the ACDM model. 
If we regard this in its conventional dark energy way, we have the dark energy with 
a constant equation of state: w = —1. If, however, we regard the cosmological 
constant as part of the gravitational sector, then we obtain 7eg*"^^^ = —QaHq/H^ 
for a matter cosmology. 

As the universe expands, it is clear from the strong energy condition, and energy 
conservation (|4.9|) , that the energy density of the universe must decrease. As the 
energy density is diluted, we ought to approach a vacuum state, which will either 
correspond to a Minkowski brane, or a self-accelerating brane. We are interested in 
the case where we approach the Minkowski vacuum. Now, in the absence of any phase 
transitions, we expect H to vary continously during the cosmological expansion. It 
follows that our brane cosmology must be consistent over a range < H < H^g,^. 

To be consistent, our cosmology must adhere to certain rules. To begin with, 
the strong energy condition requires that the energy density is non-negative. It only 
vanishes for the vacuum brane, which we have taken to be the Minkowski solution. 
We therefore require that F{H^) > over the range < H < i/max (see figure 1). It 
is easy to check that this rule implies that 

F'(0) = 6x > (4.13) 

The alternative scenario, where -F'(O) < 0, is also shown in figure ^. Here there is a 
small region close to > for which F{H^) < and therefore < 0. This is an 
unphysical regime. The only way to get a physical vacuum is for F' to change sign, 
so that we have a self- accelerating vacuum with H = Hq. The Minkowski brane in 
this instance represents an isolated vacuum which will never be approached, even as 
the energy density is diluted. 

Thus, condition ( [4.13|) guarantees that the Minkowski vacuum is the physical 



one. However, recall from the previous section that unless the radion decoupled for 
one reason or another, we required x < in order to avoid a ghost. We then arrive at 
the following important and somewhat surprising conclusion: the case of a physically 
safe, i.e. perturbatively stable, radion mode corresponds to the unphysical cosmolog- 
ical case. Therefore if we want to stick to a Minkowski vacuum then it follows that 
the radion must be absent for the solution to be consistent. This corresponds to the 
strongly coupled case (x = 0), and the cases for which decoupling occurs (x oo). 
We will study the corresponding cosmological solutions for these "radion-free" cases 
in the next section. 
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Figure 1: An example of two possible func- 
tions p = F{H^) where there is a Minkowski vac- 
uum (the dotted line), and where the Minkowski 
vacuum is isolated (solid line). These examples 
correspond to the normal and self-accelerating 
branches of pure DGP theory respectively. 



The link between a consistent 
cosmology, and the absence of a per- 
turbative ghost can be traced back 
to the absence of a tensor zero-mode. 
This means that the large distance 
or late time behaviour is dominated 
by the scalar radion (if it exists). 
Now, if one uses the Gauss- Co dazzi 
equations to find the projection of 
the Einstein tensor on the brane , 
then after linearising, we find that 

6 (7^^, - \n^^,) = ^T^, (4.14) 



like a dark source of radiation 36 



We have left out the explicit con- 
tribution from the bulk Weyl tensor 
since it is only expected to behave 
and is actually zero for the cosmological branes 



considered here. Now the effective equation of motion ( [4.14| ) would follow from an 
effective action of the form 



S, 



eff 



(4.15) 



At late times, we now have a natural interpretation for x'- it is the effective 4D 
Planck scale. Taking this to be positive, it is clear that we will have a consistent 
AD cosmology. However, in the absence of a tensor zero mode, the conformal mode 
has nothing to mix with in the far infra-red, giving rise to a dangerous "conformal 
ghost" . 

Note that the condition -F'(O) > is actually required for the whole cosmology 



F'{H^) > 0, ioTO<H <H„ 



(4.16) 



otherwise the Strong Energy Condition will fail to hold^. The condition ( [4.16[ ) also 
guarantees that we never enter a phase of super-inflation, with H > 0. This can 
easily be checked using equation (f4.7|) and energy conservation (|4.9| ). Super- inflation 
is usually associated with phantom cosmologies, which will typically lead to ghost-like 



instabilities 37 



Last, but not least, we must demand that our cosmology passes through a stan- 
dard AD phase at some earlier time, otherwise we will run into problems with nu- 



cleosynthesis |3^. If -ffmax is sufficiently large, then when H 



i^max, the induced 



^To see this, note that our requirement that F' is positive in a neighbourhood of the origin 
means that if F' becomes negative, then it must have a zero at finite H. This in turn implies that 



{p + p) = 0, from (4.9), thus violating strong energy. 
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curvature on the brane will dominate, so that p ~ Qrn^^H^. For the accelerating 
cosmologies to be discussed in the next section, we will find that -ffmax can be taken 
to be arbitrarily large. 

5. Radion-free cosmologies 

We saw previously that the only cosmological solutions that were universally consis- 
tent were those for which the radion field was absent from the vacuum perturbations. 
This included the conformal limit (x = 0), and the decoupling limits. The latter cor- 
responds to two possibilities: a generalised Randall- Sundrum scenario or a Minkowski 
bulk on at least one side of the brane. We will look at both of these cases, as each 
contains an example of an accelerating cosmology. 

5.1 Acceleration with a decoupled radion 

The generalised Randall-Sundrum scenario corresponds to the case where we retain 
the AdS horizon on either side of the brane. The scalar excitations about the vacuum 
decouple. We also have a consistent cosmology for < < oo, defined by 

= 6mjii^2 + i2(M^ (yiPT¥ - ) (5.1) 

One can easily prove that C^H"^) = F/F'H'^ > 1 for all values of H > 0. As discussed 
in the previous section, acceleration can only occur for ordinary matter satisfying the 
strong energy condition if C{H^) falls below one. We conclude that we can never 
enter a phase of cosmic acceleration in the generalised Randall-Sundrum scenario for 
ordinary forms of matter. This comes as no surprise, since the presence of a graviton 
zero-mode guarantees standard 4D behaviour at large distances, and prevents any 
infra-red modification of gravity. 

Now consider the case where we have a Minkowski bulk on at least one side of 
the brane, and without loss of generality take k2 = 0. Again, all scalar excitations 
about the vacuum decouple. Having not yet specified the values of 9i, 02, and fci, 
the cosmology is, in general, defined by 

F{H^) = Gml^H^ - 6Mf9i (^^ H"^ + kj - k^^ - 6MI92H (5.2) 

If both sides of the bulk are Minkowski, i.e., ki = 0, then 

F{H^) = GmliH^ - 12{M^e)H (5.3) 

this is the DGP model 0. Whenever {M^6) = 0, we see that the brane cosmology 
receives no contribution from the bulk, and behaves exactly as for 4 dimensional 
GR, without a cosmological constant. As is well known, this will not give us any 
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acceleration. Now consider what happens when {M^6) ^ 0. For small if > 0, we 
have 

F[E'^) ~ -l2{KPe)H + 0{H^) (5.4) 

Since p > 0, we clearly require {M^9) < 0, otherwise the fiat vacuum would represent 
an isolated vacuum, as discussed in the previous section. However, note that when 
{M^6) < 0, we have C{H'^) > 1, for all values of if > 0, and there is never cosmic 
acceleration for ordinary forms of matter. 

The first example of an accelerating cosmology occurs when we have Minlcowsl^i 
space on one side of the brane, and AdS space on the other. This corresponds to 
the case fci 7^ ^2 = 0, with the cosmology defined by (|5.2|) . Now for small H > 0, we 
have 

F{H^) ~ -mle2H + 0{H^) (5.5) 

In order to avoid isolating the fiat vacuum, we clearly require that 62 = —1- If, in 
addition, we tak;e = — 1, it is easy to check; that C{H'^) > 1 for all values of if > 0, 
and so there can never be any acceleration for ordinary forms of matter. Perhaps we 
should not be surprised that acceleration is impossible when the warp factor decays 
away from the brane on the AdS side, since this will induce a degree of localisation. 
In contrast, when the warp factor grows away from the brane, we might expect more 
interesting dynamics, since the graviton will want to localise away from the brane, 
near the AdS boundary. This does indeed happen, and, as we shall now demonstrate 
in detail, cosmic acceleration can occur for ordinary matter if we set 61 = 1. 

In the rest of this subsection, we will focus our attention on the case where 
the AdS boundary is included in A^i, and we have Minkowski space in A^2- The 
resulting cosmology is defined by 

F{H^) = 6ml^H^ - 6M^ (^-^ + kf - ki^ + 6M^H (5.6) 

Now if Ml < M2, it is relatively easy to check that C(ii^) > 1 for all ii > 0, and 
as such, acceleration can never occur for ordinary matter. We can understand this 
result as follows. Although the warp factor grows away from the brane in Mi, the 
graviton is not so strongly localised on the AdS boundary, since a small value of Mi 
makes it easier for the graviton to propagate towards the brane. Thus, the degree of 
delocalisation away from the brane on the AdS side is not as severe as it might have 
been for larger values of Mi. 

Now let us focus in detail on the case where the degree of delocalisation is more 
severe, so that cosmic acceleration is maximized, by taking Mi > M2. It follows that 
F'{H'^) has one local minimum in ii > 0, where it takes the value 

fL„ = 6™?,-3^(l--g|) (57) 
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Recall that the cosmological phase of interest corresponds to F'{H^) > for nonzero 
H. Therefore, if F^in < 0, we can only have a consistent cosmology for < if < i^max 
where if max, is finite. For < H < i^max , we can use the fact that F'^H"^) < 
to show that C is an increasing function, and so C{H'^) > C(0) = 1, thus cosmic 
acceleration can never occur for ordinary matter if F^j^ < 0. In contrast, if F^^j^ > 0, 
then it follows that F'{H'^) > for all H > 0. There is only one cosmological phase, 
so we can take i^max to be infinite. This condition on F^^^ enables us to place the 
following lower bound on 



pi 



Given that our cosmology is consistent for arbitrarily large values of H, it is instruc- 
tive to consider the asymptotic behaviour of C: 



C = l- 



Mf - M| 



H-^ + 0{H-^) (5.9) 



There are two things to note here. First, at very large if, C ~ 1 to leading order, 
so the standard AD cosmology is reproduced at early times. This is due to the 
induced curvature term on the brane dominating the UV behaviour. Second, the 
first correction demonstrates that C starts to fall below one, since Mi > M2, which 
is precisely the sort of behaviour we hope to see in order to get cosmic acceleration 
from ordinary matter. Quite how much acceleration can be obtained depends on the 
nature of the cosmological fluid (the value of w), and the minimum value of C{H'^). 
It is the latter that measures the degree to which modified gravity is contributing to 
the acceleration. 

To maximize cosmic acceleration, we therefore require the minimum possible 
value ofC{H^y. 



^ ^'^^r) (cos^.9 + e)J/^-2(Vl + g^-l) + 2g.in^ff (5.10) 

(cos3 (5 + e)ii'2 - H^/\/l + H"^ + H sin^* [3 

where we have defined 

with < /3 < 7r/2, and e > 0, which is consistent with Mi > M2 and the bound 
( |5.8|) . Unfortunately there is no simple analytic minimization of this functional, 
however, it is not difficult to establish that the minimum occurs along e = 0, i.e. on 
the boundary of our allowed range of the 4-dimensional Planck mass. Therefore, we 
can never actually attain the maximal value of acceleration, however, numerically. 
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we see that the minimum values one can typically achieve are of order C ~ 0.43. For 
a matter dominated universe {w = 0), we can therefore choose /3 and e such that the 
effective equation of state, 7eflf, falls as low as ~ —0.57, and for radiation {w = 1/3) 
~ 0.43. 

To see if this cosmological model is viable, it is not sufficient to demonstrate 
that the effective equation of state bottoms out at some reasonable negative value, 
as it may be that the actual cosmological evolution does not spend sufficient time in 
this negative region to have a significant era of late time acceleration. We therefore 
need to track the scale factor throughout time to demonstrate that the period and 
amount of acceleration is at least potentially sufficient. Differentiating the Friedmann 
equation gives the cosmological evolution of H: 

^-§ = -l{l + w)H^C{H^) (5.12) 

where t = kit. This is, in general, a somewhat involved NLDE for H, and is best 
integrated numerically. We illustrate acceleration in this model by taking a specific 
example, for which C is shown in figure |^. There are a number of generic features 
apparent in this example. At early times, when H is large {a ~ 7r/2), C is close to 
one, and so we are in the standard iD regime, because the induced curvature term 
dominates in the UV. As the universe starts to expand, and cool down, then C starts 
to fall below one. If we are in a period of matter domination {w = 0), we enter an 
accelerating phase when C falls below 2/3. The minimum value of C is around 0.5, at 
which point we have maximal acceleration, whatever the cosmological fluid. C then 
starts to increase, and we exit the accelerating phase for good. A maximum value 
of C occurs when H (3. Here we are in the phase of greatest deceleration. This 
maximum value actually diverges when e = 0, at the point at which F' vanishes. 
Beyond this phase, at very late times, C approaches its limiting value of 2, and hence 
the cosmology becomes that of a "stiff matter" universe, where p = p- 
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The evolution of the scale factor can be found by integrating (|5.12|) , and the co- 
moving Hubble radius is plotted in figure ^ for both matter and radiation dominated 
cosmologies. This shows that while the matter cosmology experiences a significant 
era of acceleration, the radiation cosmology has only a period of marginal accelera- 
tion, which would be more accurately described as 'coasting' pl| . 




10 20 30 40 ^ 10 20 30 40 



Figure 3: Plot of the comoving Hubble radius as a function of time, measured in units of 
k^^, for a) a matter cosmology, and b) a radiation cosmology, with the same values of /3 
and e as in figure ||. 

Focussing on the matter cosmology, from the plot of the comoving Hubble radius, 
we see that we enter a period of acceleration at around t ~ 0.5/cj~^, and remain in 
an accelerating phase until t ~ 32ki^. Since we chose a flat universe, our value of 
ki is determined by the current value of Qm, and the Hubble parameter Hq. Since 
we are using this model for illustrative purposes, we have not performed a best fit 
of the cosmological parameters in our model to the data, however, choosing some 
reasonable values of ki gives the luminosity-redshift plot of figure ^. 

Clearly the model gives a plausible fit to the data, however, more work needs to 
be done to establish its viability, as well as its comparative merits to other models 
[0. In particular, we have taken standard typical values of Hq and Qm, which of 
course are the result of a fit of Einstein evolution to various different data sets. It 
is possible that using the evolution predicted by our model, the best fit values of 
these parameters may be different, which was the motivation for giving a range of 
commonsense values of ki. Even so, it does seem as if our model gives a better fit 
than fiat DGP, even for the worst case scenario. 

5.2 The conformal cosmology 

The remaining "radion-free" cosmology is interesting, because it leads to the greatest 
amount of acceleration. This is the conformal case, corresponding to 

X = ml,-^-^J = (5.13) 
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Figure 4: Magnitude redshift plots of our cosmology (red curves) for a range of realistic 
ki compared to: ACDM (0^ = 0.3, Oa = 0.7) (blue), open matter (^7^ = 0.3, = 0) 
(green), and the DGP model (purple). The second plot shows a subtraction of the predicted 
magnitudes from ACDM. The labels on the red plots show the values of the parameter ki 
as a multiple of Hq. 

In this limit, an extra symmetry kicks to linear order in perturbation theory that 
eliminates the radion degree of freedom. It also prevents coupling, at linear order, 
to the trace of the energy momentum tensor on the brane. If a non-zero trace is 
present, the geometry responds non-linearly. It is precisely this non-linear effect that 
gives rise to the acceleration at late times, as we shall soon see. 

The first thing to note is that the cosmology is consistent for all values of if > 0, 
provided > and F"{0) > 0. At early times the induced curvature on the brane 
guarantees that we experience the standard AD evolution. The extra dimensions then 
begin to open up and we enter an accelerating phase at later times. To undertand 
the late time behaviour, let us approximate equation (^?7|) by a Taylor expansion 
near H = 

p = H^F'iO) + —F"{0) + — F'"(0) + O (H^) (5.14) 
2 6 

where 

F'(0)=6x, no) = 3^^^, ^0) = --^^^ (5.15) 

For X 7^ 0, p is linear in to leading order. Cosmic acceleration in this instance 
is impossible. It is only as we take the conformal limit, x ^ 0; that the linear term 
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goes away, and acceleration becomes possible. This is because p is either quadratic, 
or even cubic in at late times. This will lead to power law acceleration, since 

pocH^"" =^ a(r) ~ r^^rr^ (5.16) 

The quadratic dependence (n = 2) occurs when F"{0) > 0, or equivalently 




(5.17) 



Note that we cannot have 6i = 62 = —1, which comes as no surprise since this is just 
the generalised Randall- Sundrum scenario discussed previously, for which gravity is 
localised on the brane. To get an idea of how much acceleration occurs at late times, 
note that for small H, C{H^) = | + 0{H^) and so a pressureless cosmological fluid 
will lead to an effective equation of state, 7eff = —1/2. 

The cubic dependence (n = 3) occurs when F"{0) = 0, or equivalently 




(5.18) 



In order to avoid having an isolated vacuum at if = we must also require that 
F'"(0) > 0. This scenario leads to even more acceleration: for small H, C{H'^) = 
^ + 0{H'^) and so a pressureless fluid has an effective equation of state 7efT = —2/3 at 
very late times. Naively, we might have expected the greatest acceleration to occur 
when the AdS boundary is retained on both sides, but this turns out not to be the 
case here. If we assume, without loss of generality, that Mi > M2, then we must 
also have 61 = 1, 62 = —1. This means that in Aii, we retain the AdS boundary, 
whereas in A^2 we retain the AdS horizon. 

We can understand this puzzling result as follows. For these conformal cosmolo- 
gies, the scalar radion is strongly coupled and therefore plays a more important role 
than the spin 2 graviton mode. While the graviton likes to localise where the warp 
factor is greatest, the opposite is true for the radion. If the radion is dominant, as 
is the case here, it is clear that gravity on the brane is weakest when the radion 
localises away from the brane, near the AdS horizon. Note that in the previous case 
of a Minkowski bulk, the radion was weakly coupled, and so the graviton played the 
dominant role. 

For comparison with the Minkowski model, we choose a specific example with 
M2/M1 = ^2/^1 = sin 0.6, and show the comoving Hubble scale and a similar sub- 
tracted luminosity redshift plot in figure ^ Once again, we have not performed a 
"best-fit" analysis, but simply taken two reasonable values of the model parameters. 

In principle, if we could make p quartic in we could achieve yet more late 
time acceleration. Unfortunately, this would require us to set F'{0) = F"(0) = 



- 24 - 




Figure 5: The comoving Hubble scale for the conformal cosmology, and the subtracted 
redshift plot (relative to ACDM) for ki = Hq, 2Hq. Once again, H = H/ki, and t = kit. 

F'"(0) = 0, which corresponds here to the trivial case F{H'^) = 0. One interesting 
possibihty is that the addition of a Gauss-Bonnet term would allow this. The Gauss- 
Bonnet term is a higher order curvature correction to the gravitational action which 
is ghost-free (see for discussions in the context of braneworlds, and ^ for 



the boundary terms appropriate to Gauss-Bonnet). We can therefore view it as a 
possible first correction in a UV completion of the gravitational effective action. Such 
UV completions are sometimes proposed as a remedy to the problems of the pure 
DGP model. In this context, the effect of the GB term is to modify the relation 
between the constant /c, and the cosmological constants on each side of the brane: 



Aaikf = 1 ^ y/l + 4aiAi/3 (5.19) 

where we allow a to take different values on either side of the brane, in keeping with 
the asymmetric set-up. The brane tension is also modified to: 

a = -12{M^9k{l - ^ak'^)) , (5.20) 



The modified Friedmann equation is again (|4.7| ) but with F given by 



F{H^) = 6mliH^-12 (m'^O (yiPT¥ - A;) (1 - -ak^)^-32 (^M'^OaHWH^ + 

(5.21) 

From this, we now see that the Gauss-Bonnet term does not significantly modify 
our conclusions, although we cannot push the conformal model to stronger late time 
acceleration, since this is incompatible with good UV behaviour of the cosmology. 
Thus, we are driven to the interesting conclusion that UV corrections in this case 
appear to dampen acceleration, although we have not performed a full parameter 
search to check if there are regions of the GB theory with greater acceleration. 
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6. Discussion 



Up until now, braneworld models have offered only one new mechanism for reproduc- 
ing cosmic acceleration. These are the self-accelerating de Sitter like solutions found 
for DGP branes that have been shown to be unstable at the level of linearised 
perturbations |jl2|]. In this paper we have shown that a dynamical alternative mecha- 



nism exists whereby the vacuum brane is Minkowski, but when matter is present the 
universe can undergo a period of acceleration. Quite how much acceleration occurs 
depends on the details, but we have demonstrated two types of model for which there 
is a convincing amount of acceleration. 

The idea of ordinary matter driving accelerated expansion is not new, for exam- 



ple, the Cardassian model of Freese and Lewis, |^2[, obtained such acceleration by 
empirically modifying the Friedmann equation, sending p — » p + cp", where n < 2/3 
was some (unknown) parameter. While an interesting empirical model, Cardassian 
expansion did not provide an explanation of why the matter source contributed to 
the Friedmann equation in this way. In our model, we use concrete five dimensional 
physics to give rise to four dimensional modifications, where it is the geometry which 
is altered, rather than the matter source. Of course, one could always rewrite the 
modified Friedmann equation as 

H^ = Hp) (6.1) 

by algebraically solving p = F{H'^), however, as this involves the solution of a quartic 
in either H or H^, it is both messy, and one also has to take great care not to pick 
up fictitious additional solutions from branch choices. Nonetheless, one could view 
our models as providing a sound theoretical basis for an effective Cardassian model. 



For example, the Minkowski- AdS model illustrated in section ^]T] is approximately 
equivalent to Cardassian expansion with n = 0.5. 

The key feature of our braneworld model for modified gravity is that we have 
shown that our vacuum solutions are perturbatively stable. The requirement of 
perturbative stability is obviously crucial. Typically, instabilities appear in the scalar 
sector so we have focused our attention on scalar perturbations about the Minkowski 
brane vacuum. It turns out that there is a curious interplay between physically 
acceptable cosmologies, and stable vacua: if the cosmological solution is well behaved 
close to the vacuum, then the radion mode is a ghost (and vice- versa!) This surprising 
result can be understood as stemming from the "conformal ghost", as discussed in 
detail at the end of section |[ Perhaps more importantly, it greatly restricts the set of 
physically acceptable set-ups: we can only consider those cases for which the radion 
either decouples completely (x — oo), or is eliminated by some extra symmetry as 
in the conformal case (x = 0). 

The conformal cosmology discussed in section p.2| leads to the greatest amount 
of acceleration. Although we prefer to regard it as a stand alone special case of 
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enhanced symmetry, we can gain some understanding of its dynamics by considering 
the hmit of small x 7^ 0. The first thing to note is that the radion coupling scales 
like so we might also think of the conformal case as the strongly coupled limit. 
Because the radion is so strongly coupled it dominates the dynamics over the the 
spin 2 mode. Therefore the case of greatest acceleration occurs when one has the adS 
boundary on one side of the brane, and the adS horizon on the other. The presence 
of the adS boundary ensures that the radion does not decouple, and is able to draw 
gravity away from the brane by localising close to the adS horizon! 

The case of x small and positive may also be of phenomenological interest. The 
point is that although our linearised analysis suggests that the radion is a ghost, it is 
so strongly coupled that one cannot fully trust linearised perturbation theory except 
on the largest scales. One could therefore imagine a scenario whereby that scale is 
pushed beyond (say) the size of the observable horizon. Our cosmological solution is 
a non-linear solution so we can retain the power law acceleration in the presence of 
matter or radiation. 

Naturally, our models, while appearing to give a promising fit to the SN data, 
are not necessarily preferable to ACDM. The main problem is the extension of pa- 
rameter space - we can explain late time acceleration, but at the cost of an additional 
parameter. If we suppose that our four dimensional Planck scale is given, then our 
models have in general three free parameters^: fci. Mi and M2 (however, note that 
the conformal cosmology has only two free parameters, ki and Mi, say). This is to 
be compared to the two parameters of ACDM, VL^ and ^a- Clearly then, our models 



will in general fail on selection criteria (see e.g. [^] for a discussion of this issue 
in the current context), unless significant deviations to ACDM are found at higher 
redshift. However, our conformal cosmology, for which the parameter space is not 
higher dimensional, certainly merits further investigation. 

In this paper, we have been primarily concerned with producing a consistent 
cosmology, i.e. a model which has consistent particle physics (such as no ghosts) 
and also provides an accelerating cosmology. One question we have not considered 
is whether table-top or indeed solar system / astrophysics is consistent. The DGP 
model has received a great deal of attention, mainly because it appeared to provide 
a ghost-free way of getting the gravitational phenomenology first discovered in the 
GRS model, however, more careful analysis uncovered a number of difficulties (such 
as strong coupling) and a fundamental flaw: the existence of ghosts. Here, we have 
resolved the fundamental flaw by introducing an asymmetric braneworld, but we 
have not provided a solution (or indeed an examination) of any of the difficulties. 
For example, strong coupling, which is an issue in the DGP model, is also an issue 
in our conformal model. 



^Note that we are referring here to the parametrization of the matter or geometry, rather than 
the whole gamut of cosmological parameters. 
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Asymmetric braneworlds are an under-explored territory in the world of brane 
model building, largely because the asymmetry introduces additional complications 
at the level of the equations of motion. The main rationale for considering symmetric 
braneworlds is that these correspond to orbifold compactifications in string theory. 
However, the string landscape is also a popular concept and in the landscape we 



can expect many different vacua, with many different local cosmological constants 
[0, and indeed (through warping |^) the possibility of different local effective 
Planck masses. Therefore, although the idea of having different Planck masses and 
cosmological constants on either side of the brane may seem a little strange at first 
sight, it is in fact a natural, if not generic, occurance in the context of the string 
landscape. 

We end by re-emphasizing that we do not intend to present this as an improve- 
ment on the ACDM model. We merely wish to demonstrate what can be achieved 
consistently in a braneworld set-up. Ultimately, the concordance model is prefer- 
able because it provides the simplest fit to the data - i.e. the fit with the smallest 
number of parameters. However, the definition of 'simplest' should also include a 
measure of how naturally the model fits within an underlying fundamental theory 
of interactions. Whether it is easier to produce a massively fine tuned cosmological 
constant, or discrete vacua with differing parameters (or something else), will no 
doubt continue to the a source of lively debate for some time! 
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